Electromagnetic properties of the Beryllium-11 nucleus in Halo EFT 



D. R. Phillips 12 3 and H.-W. Hammer 2 

1 Department of Physics and Astronomy, Ohio University, Athens, Ohio 45701, USA 

2 Helmholtz-Institut fur Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universitat Bonn, D-531 15, 
Bonn, Germany 



Abstract. We compute electromagnetic properties of the Beryllium- 1 1 nucleus using an effective field theory 
that exploits the separation of scales in this halo system. We fix the parameters of the EFT from measured data on 
levels and scattering lengths in the 10 Be plus neutron system. We then obtain predictions for the B(E1) strength 
of the l/2 + to 1 /2~ transition in the "Be nucleus. We also compute the charge radius of the ground state of "Be. 
Agreement with experiment within the expected accuracy of a leading-order computation in this EFT is obtained. 
We also indicate how higher-order corrections that affect both s-wave and p-wave I0 Be-neutron interactions will 
affect our results. 



1 Introduction 

The first excitation of the Beryllium-10 nucleus is 3.4 MeV 
above the ground state, and that ground state has quantum 
numbers + . Meanwhile, the Beryllium-11 nucleus has a 
l/2 + state whose neutron separation energy is 504 keV, 
and a 1/2" state whose neutron separation energy is 184 
keV. The shallowness of these two states of n Be compared 
to the bound states of ll) Be suggests that they have signif- 
icant components in which a loosely-bound neutron orbits 
a 10 Be core. In this "one-neutron halo" picture the l/2 + 
is predominantly an s-wave bound state, while the 1/2" is 
predominantly a relative p-wave between the neutron and 
the core. In this paper we discuss efforts to use effective 
field theory (EFT) to systematically implement such a halo 
picture of the "Be nucleus. 

This halo viewpoint is reinforced by the fact that the 
scattering volume of n 10 Be scattering in the / = 1, j = 1/2 
channel has been determined to be [ 1 1 

a x =457±67fm 3 . (1) 

The corresponding length scale of order 8 fm is large com- 
pared to the natural length-scale of core-neutron interac- 
tions, which is » 2-3 fm. 

The datum ([TJ, together with the information on the 
bound-state energies in the 10 Be and "Be systems, helps 
us to estimate the expansion parameter in our Halo EFT. 
This is the binding energy of the halo nucleus, as compared 
to the energy required to excite the core, i.e. B\ /Bhi ~ 
1/6. Converting this to an estimate of the different distance 
scales involved, we infer that a majority of the probability 
density of "Be occupies a region outside the 1() Be core: 
^core/^haio ~ 0.4, which is consistent with the ratio im- 
plied by the numbers in the previous paragraph. This ratio 
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of distance scales is the formal expansion parameter for the 
EFT, and since it is not particularly small, the leading-order 
calculations presented here are only a first step. Eventually 
sub-leading orders must be computed in order to verify that 
the series is converging in the expected manner. 

Here we will apply this EFT to electromagnetic reac- 
tions in the "Be system. The B(El)(l/2 + — > 1/2") transi- 
tion has recently been measured to be 

B(El)(l/2 + ^ 1/2") = 0.105(12) e 2 fm 2 (2) 

using intermediate-energy Coulomb excitation 0. This is 
consistent with the older number 

B(El)(l/2 + ^ 1/2") = 0.116(12) e 2 fm 2 (3) 

from lifetime measurements [3 |. There are also two recent 
data sets on the Coulomb-induced breakup of the "Be nu- 
cleus [4 5 1 (see also Ref. |6|). Both experiments extracted 
the excitation function dB(El)/dE as a function of the en- 
ergy of the outgoing neutron, E. For low neutron energies 
this excitation function is affected by the final-state inter- 
action in the p-waves, and can be predicted in the halo pic- 
ture [lj. The non-energy- weighted sum rule, applied to the 
data of Ref. [4|, then gives a measurement of the radius of 
the ground state of "Be of 

<r 2 ) 1/2 = 5.7(4) fm. (4) 

This is consistent with the recent atomic-physics measure- 
ment of the charge radius J7) : 

r c = 2.463(16) fm, (5) 

when one accounts for the fact that the matter radius and 
the charge radius differ by a factor of the effective charge, 
Z e ff, of the "Be nucleus, 4/11. 
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All of these measurements can be addressed within the 
Halo EFT we will use here. In this theory the s- and p-wave 
states of the Beryllium- 11 nucleus are generated by core- 
neutron contact interactions. The theory does not get the 
interior part of the nuclear wave function correct, but, by 
construction, it reproduces the correct asymptotics of the 
wave functions of these states: 



u (r) = A exp(-y r); 
«i(r) = Aiexp(-yir) 1 + 



(6) 



for the 504 keV and 184 keV states, respectively. As such, 
it is not a method that is meant to compete with ab ini- 
tio calculations of this halo nucleus (see, e.g. B8I9I ) or of 
10 Be-n scattering IflOl . But it could prove complementary 
to such computations, since Halo EFT provides a way to 
ensure that the long-distance properties of the halo are cor- 
rectly taken care of. 

The quantities yo and j\ are determined by the neu- 
tron separation energies of the states in question. At lead- 
ing order (LO) in the expansion both Aq and A\ are fixed. 
(In the case of the p-wave this is related to the theorem 
discussed by Lee at the meeting (TTJ.) At next-to-leading 
order (NLO) both Aq and A \ receive corrections whose val- 
ues must be determined from neutron- 10 Be scattering data. 



2 Halo EFT for Beryllium-11 



We use the "Halo EFT" developed in Refs. H 1 21 131 to cal- 
culate the properties of the Beryllium- 1 1 nucleus. The de- 
grees of freedom in our Halo EFT treatment are the 1() Be 
core and the neutron. The EFT expansion in this case is an 
expansion in powers of <w/Z?hi g h- Here fihigh is, e.g. the ex- 
citation energy of states in 10 Be, and so is of order a few 
MeV, and u is the energy of the photon exciting the elec- 
tromagnetic transition of interest. 



2.1 Strong piece 

In our LO calculation we include the strong s-wave and 
p-wave interactions that lead to the shallow bound states 
in the 1 1 Be system through the incorporation of additional 
spin-zero and spin-one fields: 



c lio, -\ \c + n to, H I n 

1 2M) \ 2m J 



+irt 
j 

-go[ 



7o id, + 



2M„ 



T]\ id, + 



2M nc 
(rn'c' + cr'nc] 



n\(c Vy n) - (c 1 V; n J )n j 



(7) 



Here . . . represents additional p-wave interactions neces- 
sary to maintain Gallilean invariance, while c and n are the 



"core" and neutron fields. Hence, c is a bosonic field and 
n a fermionic one. The field <x represents the s-wave state 
and 7Tj the p-wave state. 

2.2 Dressing the s-wave state 

In order to treat the shallow s-wave state in the I0 Be-neutron 
system we adopt the counting that has been successfully 
developed to treat shallow s-wave states in the nucleon- 
nucleon system H 141 1511 61 171181 . This can be implemented 
by noting that then erne coupling is dimensionful, and tak- 
ing it to be of order /?haio- Meanwhile nc loops will have a 
typical size of order l/^haio0 an d so such a counting man- 
dates the resummation of nc loops when computing the <x 
propagator. This can be achieved through the Dyson equa- 
tion shown in Fig. [T] which leads to: 



DAP) = 



A) + mlpo - P 2 /(2M„ C )] - £o-(p) ' 



(8) 



with Eo-(p) the one-loop self-energy for the cr field. 



Fig. 1. Diagrammatic representation for the Dyson equation 
which incorporates the one-loop nc dressing of the field repre- 
senting the s-wave 10 Be-n bound state in the theory. Here and 
below the dashed line represents the field for the l0 Be core, and 
the thin solid line is the neutron. The thick grey line is the un- 
dressed a propagator, and the thick black line is the dressed cr 
propagator. 



This one-loop self-energy is calculated as: 



SAP) = - 



2n 



fi + iJ2m R \p - 



2M„ 



(9) 



when computed in power-law divergence subtraction (PDS) 
with a scale y, II15I16L Here we have introduced the re- 
duced and total masses of the neutron-core system: 



Mnc = m„ + m c \ 



m R = 



(10) 



Substituting Eq. |9]) into Eq. (|8j, we can set the pa- 
rameters go and Aq by computing the s-wave neutron-core 
scattering amplitude in the theory |7]) (see Fig.|2]i: 



t (E) = g 2 DAE,0), 



(11) 



in the two-body center-of-mass frame. This is then matched 
to the effective-range expansion in this channel: 



to(E) = 



2n 



1 



m R l/a ■ 



k 2 + ik 



(12) 



1 In a suitable regularization scheme, e.g. power-law diver- 
gence subtraction 11511 61 . this is true for both the real and imag- 
inary parts of the loops. 
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Fig. 2. The single Feynman diagram that is needed (once the 
dressing of the <x field has been computed) to obtain the s-wave 
neutron-core scattering amplitude. 



producing 

Itijq 1 1 

Do-(p) = -5-2 "i ^ + re g ular - ( 1 3 ) 

m z R g 2 1 - roro p - JL + B() 

In Eq. ( fl~3] > the position of the pole is determined by the 
binding energy Bq = y^lQm R ), and yo is the positive root 
of the equation: 



1 1 2 

— + o^oro - To = 0. 



(14) 



2.3 Dressing the p-wave state 

We proceed similarly for the p-wave state. The propagator 
for this state obeys the Dyson equation depicted in Fig. [3] 
Consequently it takes the form: 



DAp) = 



A l +J ll [p Q -p z /(2M nc )]-£ n ( P y 



(15) 



Computation of the one-loop self-energy from the perti- 
nent p-wave vertices in Eq. (|7]i gives, if PDS is employed: 



m R g\k 2 
67 



-ji + ik 



(16) 



with k = ^2m R E — |p'| = |p| for on-shell scattering. Con- 
sequently we obtain: 



DAP) 



3n 



m\g\ n + 3yi p - p 2 /(2M nc ) + B] 



Here j\ = y]2m R B\ is the solution of 
11,, 

— + ~nyt +y\ =0, 



+ regular. 
(18) 



(19) 



where a\ is the scattering volume, and r\ the p-wave "ef- 
fective range", which, in fact, has dimensions of 1/length. 



2.4 Fixing parameters 

Using the experimentally known values Bo = 504 keV, 
B\ = 184 keV, we infer y = 0.15 fin -1 , and y x = 0.09 
fm _1 , which are both of the expected size 1 //?haio- 

There are two possible countings for the p-wave scat- 
tering volume a t . Here we adopt that of Ref . [ 1 3 1 , and take 

ci\ ~ ^aio^ core ' as tn ^ s re q u i res on ly one fi ne tuning in the 
parameters of the underlying theory. (By comparison, the 
counting of Ref. IfTZl . where a\ ~ 



RLR 



halo 



, , requires two such 
and ji ~ l//?haio, we 



fine tunings.) With a\ 

can neglect the y\ term in Eq. (19i at leading order and 



deduce: 



(20) 



It follows that if we adopt the central experimental value 
from Eq. ([TJ we have r\ - -0.54 fin -1 . (Propagation of 
experimental errors to the final result is straightforward but 
will not be discussed in this contribution.) Parametrically 
r\ ~ 1 /R core . At NLO Eq. (|20| is corrected to: 



y\a\ 



(21) 




Fig. 3. Diagrammatic representation for the Dyson equation 
which incorporates the one-loop nc dressing of the field repre- 
senting the p-wave l0 Be-n bound state in the theory. Once again, 
the dashed line represents the field for the 10 Be core, and the thin 
solid line is the neutron. The thin double line is the undressed n 
propagator, and the thick double line is the dressed n propagator. 

We note that, since the self-energy loop is cubically 
divergent, both parameters, A\ and g\, are mandatory for 
renormalization at LO. This time we are interested in the 
p-wave core-neutron scattering amplitude in the center-of 
mass frame: 



fi(p',p;£) 



<7iP 



PAt(£,0) 

6n p' ■ P 

m R 1/aj - ±nk 2 + ik 3 ' 



(17) 



which reduces r\ to -0.72 fin -1 , a ~ 30% correction that 
is in line with the anticipated expansion parameter of Halo 
EFT in the "Be system. 

In the s-waves the situation is more straightforward: 
there we count Oq ~ R\u\ , and ro ~ ^ C ore- In consequence 
we can set ro = at LO, and obtain from Eq. 



14) 



1 

yo = — • 

a 



(22) 



Thus, the parameters in the strong sector at LO are n, 
yo (or equivalently qq), and yj. At NLO these are to be 
supplemented by ro. 



2.5 Including photons 

Photons are then included in the Lagrangian via minimal 
substitution: 



dn -> Dp = dp + ieQAp. 



(23) 
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The charge operator Q takes different values, depending 
on whether it is acting on a c field or an n field. Q n — 
for the neutron, and below we denote the eigenvalue of the 
operator Q for the c field as Q c — Q c c. Q c = 4 in the case 
of interest here, where the core is Beryllium- 10. 

If magnetic properties are to be discussed we must also 
consider the insertion of photon interactions with the neutron- 
core system which are gauge-invariant by themselves (e.g. 
the contribution of the neutron magnetic moment). But here 
our focus is on electric properties (and form factors) at 
LO and NLO, and so it is only necessary to consider how 
the Lagrangian |7]l is affected by the substitution ( f2"3"| >. At 
NNLO in the computation of these properties operators in- 
volving the electric field E and the fields c, n, and cr which 
are gauge invariant by themselves contribute to observ- 
ables. Thus, at that order there is at least one parameter in 
the Halo EFT description of Coulomb-induced breakup of 
the Beryllium- 1 1 ground state which cannot be fixed from 
10 Be-neutron scattering information alone. 



3 Results for observables 



Using Eq. |7]l plus minimal substitution ( 23 i we obtain a 
Lagrangian that describes interactions amongst the core, 
the neutron, the ground and excited states of the 11 Be nu- 
cleus, and photons. In this section we use this Lagrangian 
to compute the form factor of the s-wave state and the El 
transition from the s-wave to the p-wave state. 



3.1 s-wave form factor 

The s-wave form factor is computed by calculating the 
contribution to the irreducible vertex for Aqctct interactions 
shown in Fig. |4] This is the only diagram it is necessary to 
consider at leading order. After the application of wave- 
function renormalization, the irreducible vertex for the A () 
photon coupling to the cr state is equal to -ieQ c G c (\q\), 
where q is the three-momentum of the virtual photon. (Such 
an interpretation is valid provided the computation is car- 
ried out in the Breit frame, where the four-momentum of 
the virtual photon q = (0, q).) A straightforward calcula- 
tion yields: 



G e (|q|)=^arctan(^! 
fH \2yo 



with / = m/M„ c = m R /M. Note that G c (0) = 1, as it 
should. For the deuteron we have / = 1/2, and this reduces 
to the LO result of Ref. lfT9l 

The charge radius of the s-wave state can be extracted 
according to: 



G c (|q|) = 1 - \{rl)q 2 + . 
6 



(25) 



and an expansion of Eq. ( 24 1 in powers of |q| then yields 




Fig. 4. The LO contribution to the irreducible vertex for an A 
photon to couple to the field representing the 10 Be-neutron s-wave 
bound state. Note that there is no diagram for the photon to couple 
to the neutron as this order, since Q„ = 0. 



Putting in the number for y$ obtained in the previous 
section we find {r 2 ) 112 = 1.73 fm from this leading-order 
HEFT computation. This is about 30% smaller than the 
atomic physics measurement of (r 2 ) 1 ^ 2 = 2.463(16) fm. 
(The numbers from the non-energy-weighted sum rule are 
similar when corrected for /.) We expect NLO corrections 
to increase the charge radius, and a shift of the size needed 
to produce agreement with experiment is entirely consis- 
tent with the nominal ~ 0.4 expansion parameter of the 
Halo EFT in this system. At NLO a careful treatment of 
current conservation, which includes an operator associ- 
ated with gauging the term ~ cr' do<r in Eq. Q, still yields 
G C (Q) = 1, but also produces an increased charge radius, as 
long as ro > 0, cf. Ref. 120J21J. The precise size of the in- 
crease is fixed once the s-wave effective range ro is known. 

The charge radius of the p-wave state is also calculable 
in the Halo EFT. NLO corrections should be smaller there 
since its binding energy, and so its typical momentum, is 
lower. To our knowledge there is, as yet, no experimental 
determination of the charge radius of this state. 



3.2 E1 transition 



Now we consider the El transition from the s-wave state to 
the p-wave state. The irreducible vertex for this transition 
is depicted in Fig. [5] Here we compute the transition for 
a photon of arbitrary four momentum k = (co, k), and the 
sum of diagrams yields -iT™ where j is the polarization 



(24) index of the p-wave state and yU that of the photon. 





(26) 



Fig. 5. The two diagrams needed for the calculation of the irre- 
ducible vertex that governs the s-to-p-state transition, r - lsl in Halo 
EFT at leading order. 
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We first observe that both the diagrams depicted in Fig. [3] 
are divergent, but that the divergences cancel, as they should 
since gauge invariance precludes us from writing down any 
contact interaction that contributes to this observable at 
leading order. Current conservation at LO can be explic- 
itly checked, and we find that, as long as both diagrams are 
considered: 

= 0. (27) 

Note that if only the long-distance El mechanism on the 
left-hand side of Fig. [5] is considered, as was done, for ex- 
ample, in Ref. [22], then current conservation is not sat- 
isfied, and it appears that some input from short-distance 
physics is needed in order to define the prediction for this 
observable. 

For real photons we have k • e = and k • q = 0, with 
q the incoming momentum of the s-wave state. With these 
restrictions we can write the space-space components of 
the vertex function as: 



fjj = 5 ri r E + kfliTu. 



(28) 



Current conservation ( |27) then provides an alternative way 
to calculate Te, it tells us that: 



tofjo = k,r, 



jl E- 



(29) 



But, for r ; o, the diagram on the right of Fig.|5]need not 
be considered, and so 

r j0 (k)~ I d 3 r— —Yij(r)e' KT — — , (30) 

where u\ and uq are the leading-order wave functions of the 
s- and p-wave states, given by Eq. dfjk. As |k| — > Eq. ( 30 1 
reduces to: 



r j0 (k)~kjf 



dru\(r)riiQ(r), 



(31) 



an equation in which, of course, the integral is the canoni- 
cal form of the E 1 matrix element. 

Using the definition of B (El) EH we find that the B(E1) 
strength for this transition is related to the renormalized, ir- 
reducible vertex by: 



B(El) = 



i r, 



\n \ a) 



where 




at leading order. 
Hence, 



B(El) 



Z 2 - 2 



(32) 

(33) 
(34) 



An -3n 



2yi + ro 



(To +n) 2 



(35) 



is the leading-order Halo EFT result. No cutoff parameter 
is needed in order to get a finite result for B(E1): our value 



is finite without regularization. We note that the result ( |3~5j ) 
is "universal" in the sense that it applies to any El s-to-p- 
wave transition in a one-neutron halo nucleus. Once r\,j\, 
and yo are known for a given one-neutron halo the predic- 
tion ( |35] l is accurate up to corrections of order R COK /Rha\o- 
In the case of "Be we insert the numerical values from 
Sec. Eland find: 



Blo(EI) 



0.085 e 2 fm 2 . 



(36) 



This is about 20% too low compared to the value Q ob- 
tained in Ref. . 

NLO corrections from the wave-function renormaliza- 
tion factors associated with the s-wave and p-wave fields 
will serve to increase the leading-order prediction. We ex- 
pect them to be of the size necessary to bring the predic- 
tion into agreement with the datum (Ej). As was discussed 
above in the case of r c , there is a prediction for B(E1) at 
NLO if we can fix the value of ro from some other data. In 
particular, we are optimistic that we can fix the magnitude 
of ro by examining experimental results for the low-energy 
El strength function in breakup to the 10 Be-neutron chan- 
nel: dB(El)/dE. The computation of that observable is a 
straightforward extension of the s-wave bound state to p- 
wave bond state El calculation carried out here. 

We emphasize that if the s-wave effective range, ro, is 
known then the NLO result is predictive. The first contri- 
bution of physics at scale R core that cannot be fixed from 
hadronic-interaction observables does not enter the result 
for B(E1) until next-to-next-to-leading order, i.e. it is sup- 
pressed by filo/fihi ~ (Acore/flhalo) 2 - 



4 Conclusion 

This brief discussion of electromagnetic observables in the 
Beryllium- 11 system already displays the significant re- 
cent experimental activity that has been focused on this 
nucleus. Coulomb excitation has been used at a variety 
of facilities to probe El transitions, and atomic-physics 
experiments have made great strides through advances in 
trapping technology. This means that the time is ripe for a 
detailed analysis of electromagnetic properties of halo sys- 
tems. Here we have shown how EFT can provide such an 
analysis. 

The LO result in the EFT mirrors elegant analytic ap- 
proaches to halo nuclei (see, e.g. Refs. [1 22 23]). But, in 
contrast to those works, there is no regulator dependence in 
the LO result for the El strength. This is a consequence of 
current conservation in our formalism. Moreover, the cor- 
rection of <9(/? cole //?haio) to the LO prediction can be sys- 
tematically calculated using EFT techniques. We note that 
there are no spectroscopic factors in our approach. They 
represent the effect of physics at scale /? core , and the or- 
der at which such "short-distance physics" impacts elec- 
tromagnetic observables can be delineated in the EFT. In 
the Beryllium- 11 system the EFT expansion parameter is 
~ 0.4, and the LO predictions for r c and B(E1) obtained 
here are in agreement with experimental data to that ex- 
pected level of accuracy. Improving this LO prediction to 
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NLO will require us to fix ro, the effective range for s-wave 
neutron-core scattering, from experiment. This can per- 
haps be done using data on dissociation of the Beryllium- 
1 1 nucleus to the continuum. 

The Halo EFT therefore provides a complement to ab 
initio methods [8 9 10 1, which can struggle to describe El 
transitions and radii in these extended systems because of 
the widely varying core and halo scales that are present in 
the problem. In Halo EFT this wide separation of scales 
is the basis for the calculation. Input that summarizes the 
physics at scale R cole can be taken from either simulation 
or experiment, and the EFT is then used to predict the out- 
come of experiments that probe dynamics at the halo scale. 
The application of this approach to other one-neutron ha- 
los, and to two-neutron halos such as "Li, are obvious next 
steps. 
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